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Abstract 

We have recently developed some simple continuum models of static gran- 
ular media which display "fragile" behaviour: they predict that the medium 
is unable to support certain types of infinitesimal load (which we call "in- 
compatible" loads) without plastic rearrangement. We argue that a fragile 
description may be appropriate when the mechanical integrity of the medium 
arises adaptively, in response to a load, through an internal jamming process. 
We hypothesize that a network of force chains (or "granular skeleton" ) evolves 
until it can just support the applied load, at which point it comes to rest; it 
then remains intact so long as no incompatible load is applied. Our fragile 
models exhibits unusual mechanical responses involving hyperbolic equations 
for stress propagation along fixed characteristics through the material. These 
characteristics represent force chains; their arrangement expressly depends on 
the construction history. Thus, for example, we predict a large difference in 



the stress pattern beneath two conical piles of sand, one poured from a point 
source and one created by sieving. 

Granular materials are microscopically heterogenous. Despite this, it is nat- 
ural to search for continuum models that can describe their static and dynamic 
behaviour. The problem of granular statics implicitly requires knowledge of the 
construction history of the medium. At the microscopic scale, the construction 
history determines exactly where every grain is, and how it has been deformed 
from its original shape. Given this information, the microscopic forces follow 
from the local contact mechanics. But such a microscopic description of granu- 
lar materials is, in practical terms, impossible and unlikely to be a useful guide 
to their macroscopic behaviour. For example it is often assumed that if elasticity 
governs the local contact mechanics, the continuum behaviour of the assembly 
as a whole must be elastic. This may be unjustified: the physics of the granu- 
lar assembly involves additional, strongly nonlinear physics — namely, whether 
each contact is actually present or not. If, as we believe, the contact network is 
an adaptive structure that has organized itself to support the specific load ap- 
plied during construction, it may obey continuum equations quite unlike those 
of conventional elastic or elastoplastic media. 

I. INTRODUCTION 

In this paper we consider assemblies of cohesionless rough particles, whose rigidity is 
sufficient that individual particle deformations remain always small. Such assemblies are 
sometimes argued to be governed by the continuum mechanics of a Hookean elastic solid 
(perhaps with a very high modulus) . But this implicitly assumes that each granular contact 
is capable equally of supporting tensile as compressive loads. For a cohesionless medium 
this is certainly untrue: cohesionless granular assemblies are therefore not elastic. The 
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question is not one of principle, but of degree - how important is the prohibition of tensile 
forces? This is not completely clear; some would argue that it represents a negligible 
effect and that an elastic model remains basically sound, so long as the mean stresses in 
the material remain compressive everywhere. However, a fully elastic granular assembly 
would be one in which grains were, effectively, glued permanently to their neighbours on 
first contact. Because the packing is microscopically disordered, it is possible that, during 
subsequent loading a significant fraction of such contacts would become tensile, even if the 
load being applied remains everywhere compressive on average. If so, the absence of tensile 
forces is a major, even dominant, factor. 

Notice that the absence of tensile forces is a distinct physical effect from the one addressed 
by most elastoplastic continuum theories of granular media (see, e.g., Ref. 0). These are 
like elastic models, but they allow for the fact that the ratio of shear and normal forces at 
a contact cannot exceed a fixed value determined by a coefficient of static friction; this is 
usually assumed to translate to a similar condition on the bulk stress components acting 
across any plane. The resulting plasticity is similar to that arising in metals, for example, 
and not related to the prohibition of tensile forces: it applies equally for cohesive contacts. 
Of course, in applying such theories to cohesionless media one should assume that the mean 
stresses are everywhere compressive: however, as emphasized above, this constraint does not 
ensure that individual contact forces are all compressive as is actually required. 

These considerations suggest a physically very different picture of granular media, al- 
ready well developed and respected in the sphere of discrete modelling. ||^-^ In this picture, 
nonlinear physics is dominant, and the contact network of grains is always liable to reor- 
ganize as loads are applied: it is an "adaptive structure". The contact network defines 
a loadbearing "granular skeleton" shown in Fig. |l](a): this is often thought of as a network 
of "force chains", or roughly linear chains of strong particle contacts, alongside which the 
other grains play a relatively minor role in the transmission of stress, p,^] If these ideas are 
true, the continuum mechanics of the material has to be thought about afresh. Since this 
widely-accepted picture of force-chains implies a microscopically heterogeneous character of 
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the contact network in the material, it is not necessarily obvious that a continuum descrip- 
tion of it is possible at all. However, we have in recent years developed continuum models 
for granular materials which, we now believe, do capture some of the physics of force chains, 
and of their geometrical dependence on the construction history. This interpretation, which 
has evolved significantly beyond the empiricism of our early work, is developed below. 

The proposal that granular assemblies under gravity cannot properly be described by the 
ideas of conventional elastoplasticity has been opprobiously dismissed in some quarters: we 
stand accused of ignoring all that is 'long and widely known' among geotechnical engineers. 
10| However, we are not the first to put forward such a subversive proposal. Indeed workers 



such as Trollope and Harr |[T^ have long ago developed ideas of force transfer rules among 



discrete particles, not unrelated to our own approaches, which yield continuum equations 
quite unlike those of elastoplasticity. [p!^ -[T6| More recently, dynamical hypoplastic continuum 



models have been developed [T^ which, as explained by Gudehus ||T8|] describe an 'anelastic 



behaviour without [the] elastic range, flow conditions and flow rules of elastoplasticity'. Our 
own models, though not explicitly dynamic, are similarly anelastic in a specific manner that 
we describe as "fragile". 

In Section |I| below, we describe a generic "jamming" mechanism for the construction of a 
granular skeleton that, we argue, points toward fragile mechanical behaviour. This scenario 
is related, but not identical, to several other current ideas in the literature on granular 
media. [P,|6|, p!7| , p!9H 25|| These include the emergence of rigidity by successive buckling of force 



chains and the concept of mechanical percolation. [|T^] In particular there is a strong 
link between fragile media and isostatic models of granular assemblies. |l|j2^ In an isostatic 
network, the requirement of force balance at the nodes is enough to determine all the forces 
acting, so these can be calculated without reference to a strain or displacement variable. 
Isostatic networks require a mean coordination number with a specific critical value [z = 2d 
with d the dimension of space). In this sense, isostatic contact networks are "exceptional", 
and may appear remote from real granular materials. 

However, it is increasingly clear pSj^ that almost all disordered packings of frictionless 



spheres actually approach an isostatic state in the rigid particle limit. Since friction is ig- 
nored, there is still a missing link between this result and the physics of real granular media 
- a link provided by the concept of force chains, as we show below (Section [11 B| ). More 
generally, the idea that the granular skeleton could engineer itself to maintain an isostatic 
or fragile state is closely connected with the concepts of self-organized criticality (soc) ||2^ 



(see also Ref. ||2^). The concepts provide a generic mechanism whereby an overdamped 
dynamical system under external forcing can come to rest at a marginally stable (critical) 
state. In the SOC scenario, this state is characterized by hierarchical (fractal) correlations 
and large noise effects (compare Fig. |l](a)). In this article we ignore these complications and 
describe only our minimalist, noise- free models of the granular skeleton; these represent, in 
effect, regular arrays of force chains. The effect of noise on the resulting continuum equations 



is of great interest, but these require a separate discussion, which is made elsewhere. |13,14 



II. COLLOIDS, JAMMING AND FRAGILE MATTER 

A. Colloids 

We start by describing a simple model of jamming in a colloid, sheared between parallel 
plates. |TB[ This is the simplest plausible scenario in which an adaptive skeleton forms in 



response to an applied load; we believe it sheds much light on the related problem of dry 



granular media as discussed in Section |T| below. We will first use it to illustrate some general 
ideas on the relationship between jamming and fragility. 

Consider a concentrated colloidal suspension of hard particles, confined between parallel 
plates at fixed separation, to which a shear stress is applied (Fig. |l|(b) and | (a)). Above a 
certain threshold of stress, this system exhibits enters a regime of strong shear thickening; see. 



e.g., Ref. ||28|. The effect can be observed in the kitchen, by stirring a concentrated suspension 
of corn-starch with a spoon. In fact, computer simulations suggest that, at least under 
certain idealized conditions, the material will jam completely and cease to flow, no matter 
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how long the stress is maintained. In these simulations, jamming apparently occurs 
because the particles form "force chains" along the compressional direction (Fig. |l](b)). 
Even for spherical particles the lubrication films cannot prevent direct interparticle contacts; 
once these arise, an array or network of force chains can indeed support the shear stress 
indefinitely. (We ignore Brownian motion, here and below, as do the simulations; this could 
cause the jammed state to have finite lifetime.) 

To model the jammed state, we start from a simple idealization of a force chain: a linear 
string of at least three rigid particles in point contact. Crucially, this chain can only support 
loads along its own axis (Fig.^ (a)): successive contacts must be collinear, with the forces 



along the line of contacts, to prevent torques on particles within the chain. [|1^] Note that 
neither friction at the contacts, nor particle aspherity, can change this "longitudinal force" 
rule. (Particle deformability, however, does matter; see Section [III Q below.) 

As a minimal model of the jammed colloid, we therefore take an assembly of such force 
chains, characterized by a unique director (a headless unit vector) n, in a sea of "spectator" 
particles, and incompressible solvent. This is obviously oversimplified, for we ignore com- 
pletely any interactions between chains, the deflections caused by weak interactions with 
the spectator particles, and the fact that there must be some spread in the orientation of 
the force chains themselves. Nonetheless, with these assumptions, in static equilibrium and 
with no body forces acting, the pressure tensor pij (defined as Pij = —c^ij, with cTj^ the usual 
stress tensor) must obey 

Pij = P5ij + AriiUj (1) 

Here P is an isotropic fluid pressure, and A (> 0) a compressive stress carried by the force 
chains. 



B. Jamming and Fragile Matter 



Eq. (|I|) defines a material that is mechanically very unusual. It permits static equilibrium 
only so long as the applied compression is along n; while this remains true, incremental loads 
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(an increase or decrease in stresses at fixed compression axis of tlie stress tensor) can be 
accommodated reversibly, by wliat is (at tlie particle contact scale) an elastic mechanism. 
But the material is certainly not an ordinary elastic body, for if instead one tries to shear 
the sample in a slightly different direction (causing a rotation of the principal stress axes) 
static equilibrium cannot be maintained without changing the director n. Now, n describes 
the orientation of a set of force chains that pick their ways through a dense sea of spectator 
particles. Accordingly n cannot simply rotate; instead, the existing force chains must be 
abandoned and new ones created with a slightly different orientation. This entails dissipative, 
plastic, reorganization, as the particles start to move but then re-jam in a configuration that 
can support the new load. The entire contact network has to reconstruct itself to adapt to 
the new load conditions; within the model, this is true even if the compression direction is 
rotated only by an infinitesimal amount. 

Our model jammed colloid is thus an idealized example of "fragile matter": it can stati- 
cally support applied shear stresses (within some range), but only by virtue of a self-organized 
internal structure, whose mechanical properties have evolved directly in response to the load 
itself. Its incremental response can be elastic only to compatible loads; incompatible loads 
(in this case, those of a different compression axis), even if small, will cause finite, plastic 
reorganizations. The inability to elastically support some infinitesimal loads is our chosen 
technical definition of the term "fragile" . |jT6[ 



We argue that jamming may lead generically to mechanical fragility, at least in systems 
with overdamped internal dynamics. Such a system is likely to arrests as soon as it can 
support the external load; since the load is only just supported, one expects the state to be 
only marginally stable. Any incompatible perturbations then force rearrangement; this will 
leave the system in a newly jammed but (by the same argument) equally fragile state. This 



scenario is related, but not identical, to several other ideas in the literature. |T7| , p!9| -|23[] The 
link between jamming and fragility is schematically illustrated in Fig. ^. 

Now consider again the idealized jammed colloid of (Fig. ^ (a)). So far we allowed for an 
external stress field (imposed a the plates) but no body forces. What body forces can the 



system support without plastic rotation of the director? Various models are possible. One 
is to assume that Eq. (|1|) continues to apply, with P(r) and A(r) now varying in space. If 
P is a simple fluid pressure, a localized body force can be supported only if it acts along 
n. Thus (as in a bulk fluid) no static Green function exists for a general body force. (Note 
that, since Eq. (|I|) is already written as a continuum equation, such a Green function would 
describes the response to a load that is localized in space but nonetheless acts on many 
particles in some mesoscopic neighbourhood.) For example, if the particles in Fig. ^ (a) were 
to become subject to a gravitational force along y, then the existing force chains could not 
sustain this but would reorganize. Applying the longitudinal force rule, the new shape is 
easily found to be a catenary, as realized by Hooke, [0 and emphasized by Edwards. 



On the other hand, a general body force can be supported, in three dimensions, if there 
are several different orientations of force chain, possibly forming a network or "granular 
skeleton" . P-|6|,p!7| A minimal model for this is: 



Pij = Ai Uirij + A2 mirrij + A3 klj (2) 

with n, m, 1 directors along three nonparallel populations of force chains; the A's are com- 
pressive pressures acting along these. Body forces cause Ai 2,3 to vary in space. 

We can thus distinguish two levels of fragility, according to whether incompatible loads 
include localized body forces ( fett/A; fragility, e.g. Eq. (|l|)), or are limited to forces acting at the 
boundary {boundary fragilitj, e.g. Eq. (Q)). In disordered systems one should also distinguish 
between macro-fragile responses involving changes in the mean orientation of force chains, 
and the micro-fragile responses of individual contacts. We expect micro-fragility in granular 
materials (see Ref. |^), although the models discussed here, which exclude randomness, 
are only macro-fragile; in practice the distinction may become blurred. In any case, these 
various types of fragility should not be associated too strongly with minimal models such as 
Eqs. (^Q). It is clear that many granular skeletons have a complex network structure where 
many more than three directions of force chains exist. Such a network may nonetheless be 
fragile. 
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Fragility in fact requires any connected granular skeleton of force chains, obeying the 
longitudinal force rule (lfr), to have a mean coordination number z = 2d with d dimension 
of space (e.g. Fig. ^ (b) in two dimensions). This coordination number describes the skeleton, 
rather than the medium as a whole; but otherwise it is the same rule as applies for packings 
of frictionless hard spheres. These also obey the lfr - not because of force chains, but 
because there is no friction. Regular packings of frictionless spheres (which show isostatic 



mechanics) have been studied in detail recently; |]I],|T3|,^ and Moukarzel has argued that 
disordered frictionless packings of hard spheres are also generically fragile |^ (see also 
Ref. ||2^). These arguments appear to depend only on the lfr and the absence of tensile 
forces, so they should, if correct, equally apply to any granular skeleton that is made of force 
chains of three or more completely rigid particles. 



C. Fixed Principal Axis Model 

Returning to the simple model of Eq. (^, the chosen values of the three directors (two 
in d = 2) clearly should depend on how the system came to be jammed (its "construction 
history"). If it jammed in response to a constant external stress, switched on suddenly 
at some earlier time, one can argue that the history is specified purely by the stress tensor 
itself. In this case, if one director points along the major compression axis then by symmetry 
any others should lie at rightangles to it (Fig. |^ (b)). Applying a similar argument to the 
intermediate axis leads to the ansatz that all three directors lie along principal stress axes; 
this is perhaps the simplest model in = 3. One version of this argument links force chains 
with the fabric tensor, which is then taken coaxial with the stress. Q (The fabric 
tensor is the second moment of the orientational distribution function for contacts and/or 
interparticle forces.) 

With the ansatz of perpendicular directors as just described, Eq. (^ becomes a "fixed 
principle axes" (fpa) model. P,p|,p!3| Although grossly oversimplified, this leads to nontrivial 
predictions for the jammed state in the colloidal problem, such as a constant ratio of the 
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shear and normal stresses when these are varied in the jammed regime. Such constancy is 
indeed reported by Laun |^ in the regime of strong shear thickening; see Ref. |]I6[ . 



III. GRANULAR MATERIALS 

We beheve that these simple ideas on jamming and fragility in colloids are equally relevant 
to cohesionless, dry granular media constructed from hard frictional particles. For although 
the formation of dry granular aggregates under gravity is not normally described in terms of 
jamming, it is a closely related process. Indeed, the filling of silos and the motion of a piston 
in a cylinder of grains both exhibit jamming and stick-slip phenomena associated with force 
chains; see Ref. [0. And, just as in a jammed colloid, the mechanical integrity of a sandpile 
entirely disappears as soon as the load (in this case gravity) is removed. 

In the granular context, a model like Eq. (0) can be interpreted by saying that a fragile 
granular skeleton of force chains is laid down at the time when particles are first buried at a 
free surface (see Fig. |]); so long as subsequent loadings are compatible with this structure, 
the skeleton will remain intact - if not grain by grain, then at least in its average properties. 
If in addition the skeleton is rectilinear (perpendicular directors) this forces the principal 
axes to maintain forever the orientation they had close to the free surface (fpa model). 
However, we do not insist on this last property and other models, which correspond to an 
oblique family of directors in Eq. (|^), will be described below. P, p^ , p^ 



In what follows we review in more detail the nature of our fragile models and the role they 
might play within a continuum mechanical description of granular media. We will mainly be 
concerned with the standard sandpile, which we define to be a conical pile, constructed by 
pouring cohesionless hard grains from a point source onto a perfectly rough, rigid support 
as shown in Fig. |^. We assume that this construction leads to a series of shallow surface 
avalanches whereby all grains have come to rest, at the point of burial, very close to the free 



surface of the pile. (Very different conditions may apply for wedges of sand; see Ref. |]T5[.) An 



alternative history is the sieved pile, which is a cone created by sieving a series of concentric 
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discs one on top of the other. In the standard sandpile, it is well known that the vertical 
normal stress has a minimum, not a maximum, beneath the apex. ||3^j3^ A striking feature 
of our modelling approach is that it not only accounts for this "stress-dip" reasonably well, 
but predicts that it should be entirely absent for a sieved pile. This proposal is currently 
being subject to careful experimental verification. 



A. Continuum Modelling of Granular Media 

The equations of stress continuity express the fact that, in static equilibrium, the forces 
acting on a small element of material must balance. For a conical pile of sand we have, in 
d = 3 dimensions, 

drarr + O^a^r = (^{(^XX ~ (^rr)/r (3) 
drarz + dz(Tzz = g - Pl^rz/r (4) 

d^aij = (5) 

where /3 = 1. Here z, r and x are cylindrical polar coordinates, with z the downward vertical. 
We take r = as a symmetry axis, so that = o'zx = 0; is the force of gravity per unit 
volume; aij is the usual stress tensor which is symmetric in The equations for = 2 are 
obtained by setting /3 = in (|^,|) and suppressing @. These describe a wedge of constant 
cross section and infinite extent in the third dimension. 

The Coulomb law states that, at any point in a cohesionless granular medium, the shear 
force acting across any plane must be smaller in magnitude than tan (p times the compressive 
normal force. Here (p is the angle of friction, a material parameter which, in simple models, is 
equal to the angle of repose. We accept this here, while noting that (i) (p in principle depends 
on the texture (or fabric) of the medium and hence on its construction history; (ii) for a 
highly anisotropic packing, the existence of an orientation-independent is questionable; 
(iii) the identification of (p with the repose angle ignores some complications such as the 
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Bagnold hysteresis effect (which may in turn be coupled to density changes). Setting these 
to one side, we note that the Coulomb law is an inequality: therefore, when combined with 
stress continuity, it cannot lead to closed equations for the granular stresses. To close the 
system of equations, further assumptions are clearly required. One choice is to assume that 
the material is an elastic continuum wherever it does not violate the Coulomb condition. 
(This is the simplest possible type of elastoplastic model.) A second choice it to treat the 
Coulomb condition as though it were an equality. This is the basis of the so-called "rigid 
plastic" approach to granular media. We return to both of these modelling schemes after 
first describing our own approach. 

1. Constitutive Relations Among Stresses 

We view cohesionless granular matter as assembly of rigid particles held up by friction. 
The static indeterminacy of frictional forces can, we argue, then be circumvented by assuming 
the existence of some local constitutive relations (c.r.'s) among components of the stress 
tensor. [0-^ The c.r.'s among stress components are taken to encode the network of contacts 
in the granular packing geometry; they therefore depend explicitly on its construction history. 
The task is then to postulate and/or justify physically suitable c.r.'s among stresses, of which 
only one (the primary c.r.) is required for systems with two dimensional symmetry, such as a 
wedge of sand; for a three dimensional symmetric assembly (the conical sandpile) a secondary 
c.r. is also needed. 

The above nomenclature has caused confusion to some commentators on our work. In 
solid mechanics the term 'constitutive relation' normally refers to a material-dependent 
equation relating stress and strain. In fluid mechanics one has instead equations relating 
stress and (in the general case of a viscoelastic fluid) strain-rate history. Instead, our models 
of granular media entail equations relating stress components to one another, in a manner 
that we take to depend on the construction history of the material. Clearly such equations are 
intended to describe constitutive properties of the medium: they relate its state of stress to 
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other discernable features of its physical state. We see no alternative to the term 'constitutive 
relations' for such equations. 

In the simplest case, which is the fpa model one hypothesizes that, in each material 
element, the orientation of the stress ellipsoid became 'locked' into the material texture at 
the time when it last came to rest, and does not change in subsequent compatible load- 
ings. This is a bold, simplifying assumption, and it may indeed be far too simple, but it 
exemplifies the idea of having a local rule for stress propagation that depends explicitly on 
construction history. At first sight the idea of 'locking in' the principal axes seems to con- 
tradict the conception of an adaptive granular skeleton which can rearrange in response to 
small incremental loads. Remember though that this 'locking in' only applies for compatible 
loads - those which the existing skeleton can support. Any incompatible load will cause 
reorganization. We therefore require that any incompatible loads are specified in defining 
the construction history of the material. 

For the standard sandpile geometry (see Fig. ^), where the material comes to rest on a 
surface slip plane, such loads do not in fact arise after material is buried. The fpa constitutive 
hypothesis then leads to the following primary c.r. among stresses: 

arr = o"2z — 2 tan (f) azr (6) 

where is the angle of repose. Eq.(|) is algebraically specific to the case of a standard 
sandpile created from a point source by a series of avalanches along the free surface. The 
conceptual basis of the fpa model is not so narrow: indeed, we applied it already to jammed 
colloids in Section ||. More generally the fpa model is arguably the simplest possible choice 
for a history-dependent c.r. among stresses; but this does not mean it will be sensible in all 
geometries. 

A consequence of Eq. (^) for a standard sandpile, is that the major principal axis every- 
where bisects the angle between the vertical and the free surface. It should be noted that in 
cartesian coordinates, the fpa model reads: 

CTxx = 0-^^-2 sign(a;) tan a^z (7) 
13 



where x = ±r is horizontal. From Eq. (|^), the FPA constitutive relation is seen to be 
discontinuous on the central axis of the pile: the local texture of the packing has a singularity 
on the central axis which is reflected in the stress propagation rules of the model. (This is 
physically admissible since the centreline separates material which has avalanched to the 
left from material which has avalanched to the right.) The paradoxical requirement, on 
the centreline, that the principal axes are fixed simultaneously in two different directions 
has a simple resolution: the stress tensor turns out to be isotropic there. See Fig. |^. The 
constitutive singularity leads to an 'arching' effect for the standard sandpile, as previously 



put forward by Edwards and Oakeshott [|l^] and others. [|Tl| , p5 | 

The FPA model is one of a larger class of OSL (for "oriented stress linearity" ) models in 
which the primary constitutive relation (in the sandpile geometry) is, in Cartesians 

CTxx = Wzz + sign(x) (8) 

with 7], fi constants. Note that the boundary condition, that the free surface of a pile at 
its angle of repose is a slip plane, yields one equation linking r] and fi to 0; thus, for 
a sandpile geometry, the OSL scheme represents a one-parameter family of primary c.r.'s. 
The OSL models were developed to explain experimental data on the stress distribution 
beneath a standard sandpile. |]32|j33| , |36[| With a plausible choice of secondary c.r. (of which 
several were tried, with only minor differences resulting), the experimental data (Fig. ^ is 
found to support models in the OSL family with r] close, but perhaps not exactly equal, to 
unity (the FPA value). This is remarkable, in view of the radical simplicity of the assumptions 
made. 

As explained by Wittmer et ai, p,^] the OSL models, combined with stress continuity 
(Eq. ^) yield hyperbolic equations having fixed characteristic rays for stress propagation. In 
fact they are wave equations; |^J^ moreover they are essentially equivalent to Eq. (0), with 
(in general) an oblique triplet of directors (these become mutually orthogonal only in the 
case of fpa). The constitutive property that OSL models describe is that these characteristic 
rays (and not, in general, the principal axes) have orientations that are 'locked in' on burial 
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of an element, and do not change when a further compatible load is applied. As demonstrated 
already in Section y, there is every reason to identify such characteristics, in the continuum 
model, with the mean orientations of force chains in the underlying material. 

Note that unless the OSL parameter is chosen so that = 0, a constitutive singularity 
on the central mentioned above for the fpa case, remains. (The characteristics are 

asymmetrically disposed about the vertical axis, and invert discontinuously at the centreline 
X = 0.) The case /i = corresponds to one studied earlier by Bouchaud et al, [0 and of the 
OSL family it is the only candidate for describing a sieved pile, in which the construction 
history clearly cannot lead to a constitutive singularity at the axis of the pile. The resulting 
'bcc' model could be called a 'local Janssen model' in that it assumes /oca/ proportionality of 
horizontal and vertical compressive stresses - an assumption which, when applied globally to 
average stresses in a silo, was first made by Janssen. |^ The BCC model predicts a smooth 
maximum, not a dip, in the pressure beneath the apex of a pile. This is what we expect, 
therefore, in the case of a sieved pile, pj 



2. Rigid-Plastic Models 

A more traditional, but related, approach is one based on the (Mohr- Coulomb) rigid- 



plastic model. [^] To find so-called limit state solutions in this model, one postulates that 
the Coulomb condition is everywhere obeyed with equality. That is, through every point 
in the material there passes some plane across which the shear force is exactly tanc/) times 
the normal force. By assuming this, the model allows closure (modulo a sign ambiguity 
discussed below) of the equations for the stress without invocation of an elastic strain field. 

This limit-state analysis of the rigid plastic model is equivalent to assuming a 'constitutive 
relation' (sometimes called 'incipient failure everywhere' 0): 

1 



COS"' I 



(9) 



sin^ + l±2sin0 yl — (cot a^r/ cTzz)"^ 

whereas the Coulomb inequality requires only that cr^r lies between the two values (±) on 
the right. It is a simple exercise to show that for a sandpile at its repose angle, only one 
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solution of the resulting equations exists in which the sign choice is everywhere the same. This 
requires the negative root (conventionally referred to as an 'active' solution) and it shows 
a hump, not a dip, in the vertical normal stress beneath the apex. Savage JIOl, however, 



draws attention to a 'passive' solution, having a pronounced dip beneath the apex, [g^ The 
passive solution actually contains a pair of matching planes between an inner region where 
the positive root of (H) is taken, and an outer region where the negative is chosen. In fact 
(see Ref. |jl5|) there are more than one such matched solutions, corresponding to different 
types of discontinuity in the stress (or its gradient) at the matching plane and/or the pile 
centre. Moreover, there is no physical principle that limits the number of matching surfaces; 
by adding extra ones, a very wide variety of results might be achieved. 

It is interesting to compare the mathematics, and physics, of Eq. (|^) with that of the 
OSL models introduced above. The rigid-plastic model yields a local c.r. among stresses; like 
OSL the resulting equations are hyperbolic. It also exhibits fragility: because a yield plane 
passes through every material point, certain incremental loads will cause reorganization. 
Therefore, anyone who defends the rigid-plastic model as a cogent description of sandpiles 
cannot reasonably object to these same features in our own models. Conversely, we cannot 
object in principle to a model in which a Coulomb yield plane passes through every material 
point. However, we still see no reason why it should be a good model; in particular we 
cannot see how to make a link between the characteristic rays in this model (which are 
always load dependent) and the underlying geometry of the contact network in the medium. 
In contrast, this link arises naturally in the OSL framework. 

3. Elastoplastic models 

The simplest elastoplastic models assume a material in which a perfectly elastic behaviour 
at small strains is conjoined onto perfect plasticity (the Coulomb condition with equality) at 
larger ones. In such an approach to the standard sandpile, an inner elastic region connects 
onto an outer plastic one. In the inner elastic region the stresses obey the Navier equations, 
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which follow from those of Hookean elasticity by elimination of the strain variables. The 
corresponding strain field is usually not discussed, but tacitly treated as infinitesimal: the 
high modulus limit is taken. It has been argued that, for a sandpile on a rigid support fpa- 
like solutions can be found within a purely elastoplastic model, at least in two dimensions. 
^ However, since these show a cusp in the vertical stress on the centreline, they imply an 



infinitesimal displacement field incompatible with a continuous elasticity across the midline. 
IQ] On the other hand, it is possible to obtain a stress dip, in an elastoplastic model, by 
assuming that the supporting base is not rigid but subject to basal sag. This explanation 
cannot explain the data of Huntley ||33| which involves an indentable (rather than sagging) 
base. |jl5[ Moreover, it would predict a similar dip for a sieved pile, unlike or own models; 
experiments on this point are now available, and suggest that indeed no dip is seen in this 



case. 
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Objections to the elastoplastic approach to modelling sandpiles can also be raised at a 
much more fundamental level. p6|JT5[] Specifically, to make unambiguous predictions for the 
stresses in a sandpile, these models require boundary information which, at least for the 
simpler models, can be given no clear physical meaning or interpretation. We return to this 
point below. 



B. Fragile vs. Elastoplastic Descriptions 

1. Problems definining an elastic strain 

In the FPA model and its relatives, strain variables are not considered. No elastic modulus 
enters, and there is no intrinsic stress scale. The resulting predictions for a conical pile 
therefore obey what is usually called radial stress-field (rsf) scaling. Formally one has for 
the stresses at the base 

aij = ghsij{r/ch) (10) 
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where h is the pile height, c = cot a and Sij a reduced stress: a is the angle between the free 
surface and the horizontal so that for a pile at repose, a = (p. This form of RSF scaling, which 
involves only the forces at the base, |^T| might be called the 'weak form' and is distinct from 
the 'strong form' in which Eq. ( |I0| ) is obeyed also with z (an arbitrary height from the apex) 
replacing h (the overall height of the pile). Our OSL models obey both forms; only the weak 
form has been tested directly by experiment but it is well-confirmed in many systems (Smid 
and Novosad Huntley [p3| ). 



The observation of RSF scaling, to experimental accuracy, suggests that elastic effects 
need not be considered explicitly. This does not of itself rule out elastic or elastoplastic be- 
haviour which, at least in the limit of large modulus, can also yield equations for the stress 
from which the bulk strain fields, and hence also the modulus, cancel. (Note that it is tempt- 
ing, but entirely wrong, to assume that a similar cancellation occurs at the boundaries of 
the material; we return to this below.) The cancellation of bulk strain fields in elastoplastic 
models disguises a serious problem in their application to the standard sandpile and related 
geometries. IT^JTH] The difficulty is this: there is no obvious definition of strain or displace- 



ment for such a construction history. To define a physical displacement or strain field, one 
requires a reference state. In (say) a triaxial strain test (see, e.g., Ref. |^2|) an initial state 
is made by some reproducible procedure, and strains measured from there. The elastic part 
is identifiable in principle, by removing the applied stresses (maintaining an isotropic pres- 
sure) and seeing how much the sample recovers its shape. In contrast, a pile constructed by 
pouring grains onto its apex is not made by a plastic and/or elastic deformation from some 
initial reference state of the same continuous medium. The problem of the missing reference 
state occurs whenever the solidity of the body itself arises purely because of the load applied. 
Thus, for the jammed colloid considered in Section |I| above, the unloaded state is simply 
a fiuid. For the sandpile, it is grains fioating freely in space. On cannot satisfactorily define 
an elastic strain with respect to either of these reference states. 

A route to defining a strain variable does however exist, so long as one ignores the 
fact that tensile forces are prohibited. In effect, one assumes that when grains of sand arrive 



at the free surface, each one forms permanent or "glued" elastic contacts with its neighbours; 
p!6| this contact network can then, by assumption, elastically support arbitrary incremental 
loads. This is an admissable physical hypothesis, though contradictory to our own hypothesis 
of an adaptive, fragile granular skeleton. We do not yet know which hypothesis is more 
correct; the test of this lies in experiment. (It does not lie in a sociological comparison of 
how physicists and engineers approach their work, as offered by Savage. [§]) 

If the "glued pile" model is correct, then a strain variable is defined from the relative 
displacement that has occurred between adjoining particles since the moment they first were 
glued together. However, the resulting displacement field, found by integrating the strain, 
is unlikely to be single- valued. Put differently, if a glued assembly is created under gravity and 
then gravity is switched off, it will revert to a state in which there are residual elastic strains 
throughout the material, even though there is now no body force acting (Fig. |^(a)). This 
is because the particle contact network was itself created under partially-loaded conditions. 
Many elastic and elastoplastic calculations, such as all those reviewed by Savage, entirely 
ignore the problem of quenched stresses, and therefore embody an implausible "floating 
model" of a sandpile shown in Fig. |^ (b). 

Note that these effects do not become small when the limit of a large modulus is taken; 
the quenched stresses remain of order the stress that was acting during formation, and can 



take both signs (tensile as well as compressive). So, if one creates a glued pile under 
gravity and then slowly switches off the body force, tensile forces will arise long before g 
has gone to zero. In this sense, the response to gravity of a cohesionless pile is completely 
nonlinear. Correspondingly, in an unglued pile, no smooth deformation can connect the state 
of a pile created under gravity with an unloaded state of the same contact network: as the 
load is removed, such a pile will undergo large-scale reorganization. 
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2. Boundary conditions and determinacy in hyperbolic models 

Models, such as OSL, that assume local constitutive equations among stresses provide 
hyperbolic differential equations for the stress field. Accordingly, if one specifies a zero-force 
boundary condition at the free (upper) surface of a granular aggregate on a rough rigid 
support, then any perturbation arising from a small extra body force (in two dimensions, 
for simplicity) propagates along two characteristics passing through this point. In the OSL 
models these characteristics are, moreover, straight lines. Therefore the force at the base can 
be found simply by summing contributions from all the body forces as propagated along two 
characteristic rays onto the support; the sandpile problem is, within the modelling approach 
by Bouchaud et al. and Wittmer et al, mathematically well-posed. There is no need 
to consider any elastic strain field and the paradoxes concerning its definition in cohesionless 
poured sand, discussed above, do not arise. 

Note that in principle, one could have propagation also along the 'backward' characteris- 
tics (see Fig. |^ (a)). This is forbidden since these cut the free surface; any such propagation 
can only arise in the presence of a nonzero surface force, in violation of the boundary con- 
ditions. Therefore the fact that the propagation occurs only along downward characteristics 
is not related to the fact that gravity acts downward; it arises because we know already the 
forces acting at the free surface (they are zero). Suppose we had instead an inverse problem: 
a pile on a bed with some unspecified overload at the top surface, for which the forces acting 
at the base had been measured. In this case, the information from the known forces could be 
propagated along the upward characteristics to find the unknown overload. More generally, 
in OSL models, each characteristic ray will cut the surface of a (convex) patch of material 
at two points; the sum of the forces along the ray at the two ends must then be balanced 
by the longitudinal component of the body force integrated along the ray (see Fig. ^ (b)). 
These models are thus "boundary fragile" . 

In three dimensions, the mathematical structure of these models is somewhat altered, 
but the conclusions are basically unaffected. The propagation of stresses is governed by a 
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Green's function which is the response to a locahzed overload; OSL models predict that for 
(say) sand in a horizontal bed, the maximum response at the base is not directly beneath a 
localized overload but on a ring of finite radius (proportional to the depth) with this as its 
axis. [0 (This could be difficult to test cleanly because of noise effects, but there are related 
consequences for stress-stress correlations which are discussed in Ref. |]T3|,|T^.) On the other 



hand, for different geometries, such as sand in a bin, the stress propagation problem is not 
well-posed even with hyperbolic equations, unless something is known about the interaction 
between the medium and the sidewalls. But by assuming a constant ratio to shear and normal 
forces at the walls, further interesting predictions can be made, for example that the total 
weight increment measured at the base of a cylindrical silo, in response to an overload on 



the top, is a nonmonotonic function of the height of the fill. |^ These predictions represent 
clear signatures of hyperbolic stress propagation and, if confirmed experimentally, would be 
hard to explain by other means. 

3. The problem of elastic indeterminacy 

The well-posedness of the standard sandpile is not shared be models involving the elliptic 
equations for an elastic body. For such a material, the stresses throughout the body can be 
solved only if, at all points on the boundary, either the force distribution or a displacement 
field is specified. Accordingly, once the zero-stress boundary condition is applied at the 



free surface, nothing can in principle be calculated unless either the forces or the displace- 
ments at the base are known (and the former amounts to specifying in advance the solution 
of the problem). The problem does not arise from any uncertainty about what to do math- 
ematically: one should specify a displacement field at the base. Difficulties nonetheless arise 
if, as we argued above, no physical significance can be attributed to this displacement field 
for cohesionless poured sand. 

To give a specific example, consider the static equilibrium of an elastic cone of finite 
modulus, which is placed in an unstressed state (without gravity) onto a completely rough, 
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rigid surface; gravity is then switched on. This generates a pressure distribution with a 
smooth parabohc hump as in Fig. |]a. (The roughness can crudely be represented by a set 
of pins.) Starting from any initial configuration, another can be generated by pulling and 
pushing parts of the body horizontally across the base {i.e., changing the displacements 
there); if this is rough, the new state will still be pinned and will achieve a new static 
equilibrium. This will generate a stress distribution, across the supporting surface and within 
the pile, that differs from the original one. Indeed, if a large enough modulus is now taken 
(at fixed forces), this procedure allows one to generate arbitrary differences in the stress 
distribution while generating neither appreciable distortions in the shape of the cone, nor 
any forces at its free surface. This corresponds to a limit F — ^ oo, m ^ at fixed Yu where 
Y is the modulus and u the displacement field at the base. 

Analogous remarks apply to any simple elastoplastic theory of sandpiles, in which an 
elastic zone, in contact with part of the base, is attached at matching surfaces to a plastic 
zone. A natural presumption for the standard sandpile might be that Yu = (that is, the 
basal displacements vanish before the high modulus limit is taken). This is consistent with 
the "glued pile" interpretation of elastic models - one assumes that glue also firmly attaches 
grains to the support as they arrive. However, the same interpretation, as shown above, also 
requires explicit consideration of quenched stresses (see Fig. ||). Note in any case that elastic 
and elastoplastic predictions for the sandpile are indeterminate, in a rigorous mathematical 
sense, if the Y ^ oo limit is taken before the basal displacements u have been specified. 

Experiments (reviewed in detail in Gates et al. |jl5|]) report that for sandpiles on a rough 
rigid support, the forces on the base can be measured reproducibly; and, although subject 
to statistical fluctuations on the scale of several grains, do not vary too much among piles 
constructed in the same way. In contrast, for any simple elastic or elastoplastic model that 
does not include a specification of the basal displacements, there is a very large indeterminacy 
in the predicted stress distribution, even after averaging over any statistical fluctuations. An 
elastoplastic modeller who believes that the experiments measure something well-defined is 
then obliged to explain why and how the basal displacements (even if infinitesimal) are fixed 
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by the construction history. Note that basal sag P,|TO| is not a candidate for the missing 
mechanism, since it does not resolve the elastic indeterminacy in these models; the latter 
arises primarily from the roughness, rather than the rigidity, of the support. An alternative 
view is that of Evesque, |^5[ who directly confronts the issue of elastic indeterminacy and 
seemingly concludes that the experimental results themselves are and must be indeterminate; 
he argues that the external forces acting on the base of a pile are effectively chosen at will, 
rather than actually measured, by the experimentalist (see also Ref. |^6|). To what extent 
this viewpoint is based on experiment, and to what extent on an implicit presumption in 
favour of elastoplastic theory, is to us unclear. 



C. Crossover from Fragile to Elastic Regimes 

We have emphasized above the very different modelling assumptions of the fragile and 
elast(oplast)ic approaches to granular media. However, we have recently shown that hyper- 
bolic fragile behaviour can be recovered from an elastoplastic description by taking a strongly 
anisotropic limit. |T6[ Moreover, the crossover between elastic and hyperbolic behaviour, at 



least for one simple model of the granular skeleton, [|l6l is controlled by the deformability of 



the granular particles. For simplicity in this section, we restrict attention to the fpa model. 

The FPA model describes, by definition, a material in which the shear stress must vanish 
across a pair of orthogonal planes fixed in the medium - those normal to the (fixed) principal 
axes of the stress tensor. According to the Coulomb inequality (which we also assume) the 
shear stress must also be less than tan0 times the normal stress, across planes oriented in 
all other directions. Clearly this combination of requirements can be viewed as a limiting 
case of an elastoplastic model with an anisotropic yield condition: 

Wtn\ < (x„„ tan$(^) (11) 

where 6 is the angle between the plane normal n and the vertical (say) and t ■ n = 0. An 
anisotropic yield condition should arise, in principle, in any material having a nontrivial 
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fabric, arising from its construction history. The hmiting choice corresponding to the FPA 
model for a sandpile is ^{9) — for 9 = {n — 2(1)) /A (this corresponds to planes where n hes 
parallel to the major principal axis), and $(6*) = 4> otherwise. (There is no separate need to 
specify the second, orthogonal plane across which shear stresses vanish, since this is assured 
by the symmetry of the stress tensor.) By a similar argument, all other OSL models can also be 
cast in terms of an anisotropic yield condition, of the form 1(7^ — Cnn tan ^(6*) | < ann tan $(6*) 
where ^{9) vanishes, and ^{9) is finite for two values of 9. (This fixes a nonzero ratio of 
shear and normal stresses across certain special planes.) 

At this purely phenomenological level there is no difficulty in connecting hyperbolic 
models smoothly onto highly anisotropic elastoplastic descriptions. Specifically, consider a 
medium having an orientation-dependent friction angle ^{9) that does not actually vanish, 
but is instead very small (< e, say) in a narrow range of angles (say of order e) around 
9 — {tt — 20)/4, and approaches (f) elsewhere. (One interesting way to achieve the required 
yield anisotropy is to have a strong anisotropy in the elastic response, and then impose a 
uniform yield condition to the strains, rather than stresses.) 

Such a material will have, in principle, mixed elliptic/hyperbolic equations of the usual 
elastoplastic type. The resulting elastic and plastic regions must nonetheless arrange them- 
selves so as to obey the FPA model to within terms that vanish as e — > 0. If e is small 
but finite, then for this elastoplastic model the results will depend on the basal boundary 
condition, but only through these higher order corrections to the leading (fpa) result. Thus, 
although elastoplastic models do suffer from elastic indeterminacy (they require a basal dis- 
placement field to be specified), the extent of the infiuence of the boundary condition on 
the solution depends on the model chosen. Strong enough (fabric-dependent) anisotropy, in 
an elastoplastic description, might so constrain the solution that it is primarily the granular 
fabric (hence the construction history) and only minimally the boundary conditions which 
actually determine the stresses in the body. For models such as that given above there is 
a well-defined limit where the indeterminacy is entirely lifted, hyperbolic equations are re- 
covered, and it is quite proper to talk of local stress propagation 'rules' determined by the 
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construction history of the materiaL Our continuum modeUing framework is based precisely 
on these assumptions. 

The crossover just outhned can also be understood directly in terms of the microme- 
chanics of force chains, at least within the simplified picture developed in Section We 
consider a regular lattice of force chains (see Fig. |^ (b)), for simplicity rectangular (the FPA 
case), which is fragile if the chains can support only longitudinal forces. As mentioned in 



Section [III Q this is true so long as such paths consist of linear chains of rigid particles, 
meeting at frictional point contacts: the forces on all particles within each chain must then 
be colinear, to avoid torques. This imposes the (fpa) requirement that there are no shear 
forces across a pair of orthogonal planes normal to the force chains themselves. Suppose now 



a small degree of particle deformability is introduced. [|T6| This relaxes slightly the collinear- 
ity requirement, but only because the point contacts are now flattened (see Fig. | (b)). The 
ratio e of the maximum transverse load to the normal one will therefore vanish with (some 
power of) the mean compression. This yield criterion applies only across two special planes; 
failure across others is governed by some smooth yield requirement (such as the ordinary 
Coulomb condition: the ratio of the principal stresses lies between given limits). The granu- 
lar skeleton just described, which was fragile in the limit of rigid grains, is now governed by a 
strongly anisotropic elastoplastic yield criterion of precisely the kind described above. This 
indicates how a packing of frictional, deformable rough particles, displaying broadly con- 
ventional elastoplastic features when the deformability is significant, can approach a fragile 
limit when the limit of a large modulus is taken. (It does not prove that all packings become 
fragile in this limit.) Conversely it shows how a packing that is basically fragile in its response 
to a graviational load could nonetheless support very small incremental deformations, such 
as sound waves, by an elastic mechanism. 

The question of whether sandpiles are better described as fragile, or as ordinarily elasto- 
plastic, remains open experimentally. To some extent it may depend on the question being 
asked. However, we have argued, on various grounds, that in calculating the stresses in a 
pile under gravity a fragile description may lie closer to the true physics. 
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IV. CONCLUSIONS 



The jammed state of colloids, if it indeed exists in the laboratory, has not yet been fully 
elucidated by experiment. It is interesting that even very simple models, such as Eq. (|l]), can 
lead to nontrivial and testable predictions (such as the constancy of certain measured stress 
ratios). Such models suggest an appealing conceptual link between jamming, force chains, 
and fragile matter. |T3,TB| However, further experiments are needed to establish the degree 
to which they are useful in describing real colloids. 

For granular media, the existence of tenuous force-chain skeletons is clear; P-|S|, p!7|P^ 
the question is whether such skeletons are fragile. Several theoretical arguments have been 
given, above and elsewhere, to suggest that this may be the case, at least in the limit of 
rigid particles. Moreover, simulations show strong rearrangement under small changes of 
compression axis; the skeleton is indeed "self-organized". Experiments also suggest 



cascades of rearrangement p7| , |3T| in response to small disturbances. These findings are 
consistent with the fragile picture. 

The standard sandpile (a conical pile formed by pouring onto a rough rigid support) has 
played a central role in our discussions. From the perspective of geotechnical engineering, the 
problem of calculating stresses in the humble sandpile may appear to be of only of marginal 
importance. The physicist's view is different: the sandpile is important, because it is one of 
the simplest problems in granular mechanics imaginable. It therefore provides a test-bed for 
existing models and, if these show shortcomings, may suggest ideas for improved physical 
theories of granular media. 

Given the present state of the data, a conventional elastoplastic interpretation of the 
experimental results for sandpiles may remain tenable; more experiments are urgently re- 
quired. ^5 In the mean time, a desire to keep using tried-and-tested modelling strategies 
until these are demonstrably proven ineffective is quite understandable. We find it harder 
to accept the suggestion [|1^] that anyone who questions the general validity of traditional 
elastoplastic thinking is somehow uneducated. 
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In summary, we have discussed a new class of models for stress propagation in granu- 
lar matter. These models assume local propagation rules for stresses which depend on the 
construction history of the material and which lead to hyperbolic differential equations for 
the stresses. As such, their physical basis is substantially different from that of conventional 
elastoplastic theory. Our approach predicts 'fragile' behaviour, in which stresses are sup- 
ported by a granular skeleton of force chains that respond by finite internal rearrangement 
to certain types of infinitesimal load. Obviously, such models of granular matter might be 
incomplete in various ways. Specifically we have discussed a possible crossover to elastic 
behaviour at very small incremental loads, and to conventional elastoplasticity at very high 
mean stresses (when significant particle deformations arise). However, we believe that our 
approach, by capturing at the continuum level at least some of the physics of force chains, 
may offer important insights that lie beyond the scope of previous continuum modelling 
strategies. 
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FIGURES 





FIG. 1. (a) Granular skeleton in two-dimensional frictional hard spheres by F. Radjai et al. 
(Figure courtesy of F. Radjai.) (b) The jamming transition in a sheared colloid. The data are from 
a computer simulation of a hard sphere colloidal suspension aX (p = 0.54 which has been strained 
to 7 = 0.22. Shown in the figure are only those spheres which have come into very close contact 
(< 10~^ radius) with at least one neighbour. As seen from the figure, the contact geometry is 
strongly anisotropic and suggests the formation of "force chains" running top left to bottom right. 
(The simulation is by J. Melrose, Cavendish Laboratory; the figure is courtesy of him.) 
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FIG. 2. (a) A jammed colloid (schematic). Black: force chains; grey: other force-bearing par- 
ticles; white: spectators, (b) Idealized rectangular network of force chains. 
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FIG. 3. (a) A force chain of hard particles (any shape) can statically support only longitudinal 
compression. Note that neither friction at the contacts, nor particle aspherity, can change this 
"longitudinal force" rule, (b) Finite deformability allows small transverse loads to arise. 
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FIG. 4. Schematic "phase diagram" of a jamming system. If, as load or density is increased 
(dashed arrow), the granular skeleton arrests on first being able to support the applied load, it 
can remain indefinitely on the "marginal manifold" separating conventional solids from liquidlike 
packings. Incompatible loads will move the system to another point on the same manifold. 
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FIG. 5. Definition of normal construction history of a pile. The grains fall down from the point 
source on the pile and roll down the slopes, which are at the repose angle ^. The height of this 
pile is H{t). Sketch of the geometry of the FPA model: The stress ellipsoid has fixed inclination 
angle its ellipticity varies from zero at the centre of the pile to a maximum in the outer region. 
The outward and inward stress propagation characteristics are indicated by short-dashed and 
long-dashed lines; these are at rightangles and coincident with the principal axed of the stress 
ellipsoid. 
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FIG. 6. Comparison of FPA model Using a uniaxial secondary closure js],^ with scaled exper- 



imental data of Smid & Novosad [32| and (*) that of Brockbank et al. which was averaged 
over three piles. Upper and lower curves denote normal and shear stresses. The data is used to 
calculate the total weight of the pile which is then used as a scale factor for stresses. The horizontal 
coordinate S = rian{(j))/H is scaled by the pile height H. 
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FIG. 7. (a) The response to a localized force is found by resolving it along characteristics 
through the point of application, propagating along those which do not cut a surface on which 
the relevant force component is specified. For a pile under gravity, propagation is only along the 
downward rays, (b) Admissible boundary conditions cannot specify separately the force component 
at both ends of the same characteristic. If these forces are unbalanced (after allowing for body 
forces), static equilibrium is lost. 
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FIG. 8. Starting from an elastic cone or wedge on a rough support, any initial stress distribution 
can be converted to another by displacements with respect to the rough 'pinning' surface (a) (b) . 
Taking the limit of a high modulus (b) (c) at fixed surface forces, an arbitrary stress field remains, 
while recovering the initial shape of the cone and satisfying the free surface boundary conditions. 
This shows the physical character of 'elastic indeterminacy' for an elastic or elastoplastic body on 
a rough support. 
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FIG. 9. (a) Quenched stresses in an elastic sandpile. Layers are added in a state of zero stress to 
a pre-existing, gravitationally loaded pile. Such a pile (if gravity is removed) will spring into a new 
shape, characterized by a nonzero internal stress field which includes tensile streses. These require 
rearrangements if the grains are cohesionless; the response to gravity is intrinsically nonlinear, 
(b) The 'floating' model of a sandpile. An unstrained, isotropic elastoplastic cone is brought into 
contact with a rough surface and gravity then switched on. This is the only construction history 
we can think of that completely avoids quenched stresses in the formation of the pile. 
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